DMRG studies on linear-exchange quantum spin models in one dimension 
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We study a class of spin-1/2 quantum antiferromagnetic chains using DMRG technique. The 
exchange interaction in these models decreases linearly as a function of the separation between the 
spins, Jij = R — \i — j\ for \i — j\ < R. For the separations beyond R, the interaction is zero. 
The range parameter R takes positive integer values. The models corresponding to all the odd 
values of R are known to have the same exact doubly degenerate dimer ground state as for the 
Majumdar-Ghosh (MG) model. In fact, R = 3 is the MG model. For even R, the exact ground 
state is not known in general, except for R — 2 (the Bethe ansatz solvable Heisenberg chain) and 
in the asymptotic limit of R where the two MG dimer states again emerge as the exact ground 
state. In the present work, we numerically investigate the even-i? models whose ground state is not 
known analytically. In particular, for R = 4, 6 and 8, we have computed a number of ground state 
properties. We find that, unlike R — 2, the higher even-i?, models are spin-gapped, and show strong 
dimer-dimer correlations of the MG type. Moreover, the spin-spin correlations decay very rapidly, 
albeit showing weak periodic revivals. 

PACS numbers: 75.10.Jm, 75.40.Mg, 75.10.Kt 



I. INTRODUCTION: 



The quantum spin chains have been a subject of intense 
research effort for many decades. 1-3 While the quan- 
tum fluctuations manifest strongly in one dimension, the 
frustration also brings out interesting possibilities such 
as the spontaneous dimerization as found for the first 
time in the ground state of the Majumdar-Ghosh (MG) 
model. 4 Apart from the theoretical interests initiated by 
the Bethe's solution of the nearest-neighbor Heisenberg 
chain, or by the exact dimer ground state of the MG 
model, there exist a variety of (quasi) one dimensional 
materials which provide strong reasons for studying one- 
dimensional spin problems. 6-8 Although quite a few of 
the one-dimensional quantum spin models are known to 
admit full or partial exact analytical solution, in general, 
one still depends upon numerical approaches to inves- 
tigate these problems. The exact numerical diagonal- 
ization has been the most standard method, but it is 
severely limited by the size. 9,10 In the wake of this, the 
density matrix renormalization group (DMRG) method 
has emerged as a valuable numerical tool of computation 
for accessing large system sizes. 11-13 Although approxi- 
mate, the DMRG is an excellent technique to study the 
ground state and the lowest few excited states of one di- 
mensional model Hamiltonians. 

The present work is aimed at studying, using DMRG, a 
class of one-dimensional quantum spin-1/2 antiferromag- 
netic models in which the exchange interaction varies lin- 
early with the spin-spin distance upto a distance of R— 1, 
beyond which it is zero. The Hamiltonian for this class 



of 'linear-exchange' models is written as follows. 

R-l 
i n—1 

Here, i is the lattice-site label, and R takes positive in- 
teger values > 2, and generates different models. The 
construction of H R in Refs. 18 and 19 was inspired by 
the structure underlying the MG model. Interestingly, 
the first of these models, H2 for R = 2, is the famous 
nearest-neighbor Heisenberg problem solvable by Bethe 
ansatz. The next model H3, for R = 3, is in fact the 
MG model whose exact doubly degenerate ground state 
is given by the following singlet-dimer configurations. 

\MGl) = [1, 2) <g> [3, 4) <g) • • • ® [L - 1, L) (2a) 
\MG2) = [2, 3) ® [4, 5) <g> • • • ® [L, 1) (2b) 

Here, = (| tiij) ~ I Utj)) /V% denotes the singlet 

state formed by the spins at sites i and j. In writing 
the two MG states above, it is assumed that we have an 
even number of total spins, L, and a lattice with periodic 
boundary condition. 

Interestingly, the doubly degenerate dimer ground 
state of the MG model, H 3 , is also known to be the ex- 
act ground state of H R for any odd value of R. Thus, 
all the odd-i? linear-exchange models are exactly solvable 
for the ground state. 18,20 It is also known that the ele- 
mentary excitations, which can be viewed as dispersing 
spin-1/2 domain walls (spinons), 21 are gapped. 19 ' 22 The 
situation for the even- R models is different, however. The 
H 2 happens to be the famous integrable model solvable 
by Bethe ansatz. It is known to have algebraic spin-spin 
correlations in the ground state and gapless excitations. 
Apart from H2 , the only other rigorously understood case 
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(at least for the ground state) is the asymptotic limit of 
i?, defined as i? ~ (D(L/2) and L — > oo (1/ being the 
total number of spins on the chain). In this asymptotic 
limit, the doubly degenerate MG ground state has been 
shown to become the exact ground state of the even-i? 
models. 19 Except for H 2 and the asymptotic limit, cur- 
rently, there is no rigorous understanding of the even-i? 
linear-exchange models. 

The main objective of the present work is to investi- 
gate the even-i? linear exchange models for finite i? > 4. 
We have used DMRG technique to numerically study 
the ground state properties of these models, in partic- 
ular the H 4 , H 6 and H 8 . We find that iii?>4 behave 
qualitatively differently from H 2 , quite unlike the odd-i? 
models, all of which behave like H 3 . For example, the el- 
ementary excitations in H 2 are known to be gapless, but 
we find the even-i? ^ 2 models to be spin-gapped. Cor- 
respondingly, the spin-spin correlation in their ground 
states decays very rapidly. Moreover, they exhibit the 
MG type dimer order in the ground state which strongly 
tends, with increasing i?, towards the asymptotic behav- 
ior. Thus, we find the even-i? ^ 2 models standing closer 
to the odd-i? models, while the H 2 stands alone. This 
is contrary to a naive expectation that all the even-i? 
models might behave qualitatively similarly, as the odd- 
i? models do. Although the spin-spin correlations decay 
rapidly for the even-i? ^ 2 models, they exhibit weak pe- 
riodic revivals. We analyze this in some detail, and find 
— 1) prominent wavenumbers, corresponding to the 
correlation-revival, pretty close to the integer multiples 
of 2tt / R (in addition to 7T, which is the dominant antifer- 
romagnetic wavenumber anyway), as expected from an 
Ising type analysis. 19 

This paper is organized as follows. In Sec. II, we briefly 
discuss the DMRG algorithm and the computation of 
various physical quantities within DMRG. In Sec. Ill, we 
present the results and discussion of the ground state 
properties of H4, H 6 and ii 8 . Finally, we conclude with 
a summary. 



II. QUICK OVERVIEW OF DMRG 

The DMRG is a very useful method which allows us 
to calculate a lot of useful quantities, such as the ground 
state energy, the energy gap and the spin-spin correla- 
tions, for large system sizes. Operationally speaking, 
it is an efficient numerical scheme for truncating the 
Hilbert space of low dimensional (especially one dimen- 
sional) quantum systems. The basic idea is to optimize 
the basis states by keeping a bunch of highly probable 
eigenstates of the density matrix at successive stages of 
renormalization. 13 

The implementation of DMRG uses a superblock struc- 
ture consisting of a system and an environment. For 
the system states denoted as and the environment 
states denoted as \j), the superblock ground state, \ip), 
can be written as: |?/>) = ■ where i and j are 



summed over the complete states of the system and envi- 
ronment, respectively. The ground state density matrix 
of the superblock is defined as p — \ip)(ifj\. The re- 
duced density matrix, p, of the system block is obtained 
by tracing p over the environment. It can be written 
as: pa> 

= Ylj^ii^t'j- Here, p u > is the matrix element 
of the reduced density matrix of the system. The 
p play a crucial role in DMRG, because a selection of 
its large eigenvalue states forms the truncated basis for 
constructing a renormalized superblock. 

There are two different DMRG algorithms in use, the 
infinite and the finite system algorithms, differing by the 
choice of the environment block in relation to the system 
block. 12 In the infinite system algorithm, the environ- 
ment is the mirror image of the system. The starting 
superblock consists of four (or more) bare sites, of which 
the two on the left-hand-side may be called the system 
and the remaining two sites define the environment. Both 
the blocks grow by one site at a time, increasing the su- 
perblock size by two sites at each iteration. Since, in 
principle, there is no limit on this growth, it is called the 
infinite system DMRG. The iterations are repeated until 
the desired number of sites is reached. 

The most important step in each iteration is to trun- 
cate the system basis by keeping only a certain number of 
high probability eigenstates of p. The system block in the 
truncated basis acts as an updated left-most site, and its 
mirror image as the right-most site. Adding two new bare 
sites between them forms a renormalized superblock and 
completes one iteration of DMRG. It is important that 
the truncated basis dimension is not too small to give 
poor accuracy. If we choose m highest eigenstates, then 
a measure of the truncation error is the weight of the 
discarded states, p m = 1 — J^iLi^i, where w^'s are the 
eigenvalues of p. Naturally, keeping more states is always 
better for accuracy. But unfortunately m can not be very 
large either, because the superblock's Hilbert space grows 
exponentially which makes the computation very hard. 

In the finite system algorithm, first one uses the in- 
finite system algorithm to build up a superblock of the 
desired size L. After that, at each RG step, the system is 
made to grow at the expense of the environment (or vice 
versa), keeping the total length L fixed. Now, the system 
and the environment blocks are asymmetrical in general. 
One needs to do a few sweeps (typically three or four) 
of growing and diminishing system size to obtain more 
accurate results (specially when the system is gapless). 

For the DMRG calculations on the Hr models, where 
the interactions go beyond the nearest neighbors, we con- 
struct the superblocks in such a manner that the two 
renormalized blocks on the left and right do not directly 
interact with each other. This is done to avoid the den- 
sity matrix errors. Hence, we have taken four bare sites 
between the left and right blocks for i? = 3, 4, 5, six bare 
sites in the middle for i? = 6,7, and eight bare sites 
for i? = 8,9. We have used open boundary condition 
throughout our calculations. We have also implemented 
the total S z symmetry, and solved H± with m = 200, Hq 
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with m = 150, and Hg, with m = 95. We find that, in 
the total S z = sector for i? = 4, the truncation error 
is of the order of 1CT 13 . For iJ 6 , it is ~ 10~ 10 , and for 
Hg it is around 10~ 8 . In total S z = 1 sector, for H4 it 
is - KT 10 ; for R = 6, it is ~ 10~ 8 ; and, for R = 8 it is 
of the order of 10~ 6 . For the odd-i? models, the typical 
truncation errors are of the order of 10 -15 in the total 
S z = sector. In total S z — 1 sector, for R = 3, it is 
- 1CT 13 ; for R = 5, it is ~ 1CT 10 ; for R = 7, it is ~ 1CT 9 ; 
and, for i? = 9 it is of the order of 10~ 7 . 



III. RESULTS AND DISCUSSION 




We now present the results of DMRG calculations for 
the linear-exchange spin models. In the calculations for 
a given Hr, the nearest neighbor interaction, R — 1, is 
set as the unit of energy. To characterize the even-i? 
models, we have computed a number of useful quantities. 
Below we present them one by one, and discuss their 
implications for the nature of these models. We have 
performed calculations on the models upto R = 8. Going 
beyond it becomes computationally difficult. In any case, 
this is sufficient to get a very clear picture of the even- 
i? 7^ 2 linear-exchange models. 



A. Ground state energy 

We benchmark our DMRG codes by computing the 
exactly known ground state energies (per spin), e g , of the 
J?2 and the odd-i? models. Of course, we get the correct 
values of \ - log2 = -0.443147 for H 2 , and -3/8 for the 
odd-i? models. The e g for R = 2, 4, 6 and 8 are presented 
in Fig. 1, together with -3/8 for the odd-i?. It is clear 
that the e g for i? ^ 2 rapidly tends to -3/8. In fact, for 
i? = 8, it is already very close. The data is in agreement 
with the exact diagonalization results on small clusters, 
and clearly approaches the exact asymptotic limit. 
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FIG. 1. Ground state energy per site vs. R for L — 300. 



FIG. 2. Spin-gap vs. l/L up to L = 300. 



B. Spin-gap 

The odd-i? linear-exchange models are known to be 
spin-gapped. Also known is the gapless behavior of H^. 
The challenge is to find out if the even-i? ^ 2 models arc 
gapped or gapless. It is not obvious a-priori, although the 
behavior of e g in Fig. 1 suggests that, for large enough i?, 
the even-i? models may show up gapped behavior. To our 
surprise, we find a robust non-zero energy gap for an i? 
as small as 4, as shown in Fig. 2. The gap increases with 
increasing R. The thermodynamic values of the spin- 
gaps are estimated to be A^ =4 w 0.17, A^, =6 w 0.35 and 
A^ =s s» 0.45. The energy gap to spin excitations is esti- 
mated by computing = E^ z=1 — Eg* =0 for different 
values of L, and extrapolating to l/L = 0. Here, £^ z=0 
and Eg- =1 are the lowest energy eigenvalues in the total 
S z = and 1 sectors, respectively. 

According to the Lieb-Schultz-Mattis (LSM) theorem 
for spin-1/2 chains, the existence of an energy gap in the 
thermodynamic limit (L — > 00) implies a (doubly) de- 
generate ground state which is typically facilitated by 
an enlarged (dimerized) unit cell. The odd-i? linear- 
exchange models clearly satisfy the LSM theorem. So 
does H 2 . The LSM actually proved for H 2 a unique 
ground state and gapless excitations. 2 ' 5 In view of the 
LSM theorem, the existence of spin-gap in the even- 
i? 7^ 2 models strongly suggests dimerization in the ther- 
modynamic limit, even for finite R. It doesn't mean 
that the ground wavefunctions for these models will be 
same as the exact MG states. Nonetheless, it is expected 
that they show similar singlet-dimer order as in the MG 
state, which in the asymptotic limit is known to become 
exact. 19 



C. Dimer order 

Guided by the LSM theorem, and the asymptotic limit, 
we investigate dimerization in the ground states of iJ 4 , 
Hq and Hg,. The results of DMRG calculations for the 
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FIG. 3. Dimer-dimer correlation for L = 40. 



dimer-dimer correlation, D(r) — (Si • S2 Si +r • S2+ r ), 
are shown in Fig. 3. The data computed for L = 40, 60 
and 80 essentially looks the same. The D(r) measures 
the correlation between the total-spin states on different 
nearest-neighbor bonds. If two such bonds are perfect 
singlets, then D is equal to (-3/4) 2 =9/16 (=0.5625). 
For the MG states, D(r) switches between and 9/16 
every time r increases by one, because every alternate 
pair of nearest-neighbor spins forms a singlet. Notably, 
for R = 4, 6 and 8 too, the same qualitative behavior is 
found. Quantitatively, for H4, the minimum of D(r) is 
quite close to zero and the maximum is only slightly less 
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FIG. 4. The spin-spin correlation, C(r), vs r. 
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FIG. 5. G r = i log[(-l) r C*(r)] vs. r for L = 40. 



than 9/16. For R = 6, the numbers get pretty close to the 
ideal values. For H$, the D(r) is almost exactly like in 
the MG ground state. This data clearly demonstrates the 
nearest-neighbor singlet formation in the ground states 
of the even-i? 2 models, and puts them qualitatively 
alongside the odd-i? models. 

It should be noted that on an open boundary chain, for 
the odd-i? models, of the two degenerate MG states in 
Eqs. (2a) and (2b), only \MG1) is selected as the ground 
state. So, while the nature of the ground state remains 
unchanged, the degeneracy is 'artifically' changed to one. 
The same is happening for the even-i? ^ 2 models. We 
are getting the MG dimer order in the ground state, but 
because of the open boundary condition of the DMRG 
computations, we do not get two-fold degeneracy. 



D. Spin-spin correlation 

We have also computed spin-spin correlation, C(r) — 
(Si • Sj-j-f), which is an important characterizer of a spin 
system. For better accuracy of C(r), we have used finite 
size DMRG. The calculations have been performed on 
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FIG. 6. (— l) r C(r)e r ^ and its Fourier transform (inset). 
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FIG. 7. The structure factor 5(fe) vs fc. 



spin chains up to 100 sites long. The data for C(r) is 
presented in Fig. 4. 

The notable features of the spin-correlation data for 
even-i? ^ 2 are the following. Unlike for H2, the C(r) 
decays very rapidly to zero upon increasing r. It prompts 
us to empirically introduce a correlation length £ such 
that C(r) ~ e~ r /'>, which is consistent with the gapped 
behavior of the models. The sign of C(r) follows (— l) r , 
as generically expected for an antiferromagnet. We have 
estimated £ by computing G r = =^ log [(— l) r C(r)] such 
that Gqo = l/£. Fig. 5 clearly shows G r converging to a 
finite value for large r. The estimated numerical values 



of £ for R = 4, 6 and 8 are 4.0988, 2.82876 and 2.82876, 
respectively. One can also estimate £ from the spin-gap 
as £ ss 1/A. 24 Through this approximation, we find £ w 
5.84604 for i? = 4, £ 2.83135 for i? = 6 and £ « 2.21435 
for R = 8. Both the estimates give comparable values of 
£. The accuracy for larger i? is lesser however. 



1. Correlation-revival 

A careful look at C(r) shows weak but repeated re- 
vivals of C(r), around r ~ R, 2R and so on, of an other- 
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FIG. 8. The structure factors of the H2 and H3 models. 



wise rapidly decaying correlation function. This property 
of the spin correlations has been argued to be a signa- 
ture of the corresponding Ising problem which has an 
order with enlarged unit of size i? (a generalized sort of 
Neel state for an even-i? Ising model). 19 Here, we analyze 
this feature in some detail. 

To clearly see the correlation-revival features, we undo 
the rapid decay of C (r) by multiplying it with this empir- 
ical factor of e~ r /'% using the estimated values of £. The 
plots of e~ r ^C(r) vs. r, in Fig. 6, nicely show the repeti- 
tive feature. Since i? is roughly the period of correlation- 
revival, we expect k ~ 2 -j=? (integer n < R/2) to show up 
as special fc-points in the corresponding structure factor, 
defined below. 25 

S(*) = ~ + £££Q(r)cos*T (3) 

2=1 r=l 

Here, C;(r) = (S; • S/ +r ) is the correlation between the 
spins Si and S; +r . The computed structure factors are 
shown in Fig. 7. For any R, the S(k) always peaks at 
k = 7T, due to the factor (—1)'* in C(r). In contrast to 
the smooth S(k) for H 2 and H 3 (see Fig. 8), there exist 
identifiable shoulder- like features for H4, Hq and Hg at 
/c-points other than ir. 

As shown in Fig. 7, the S(k) for R = 4 has a clear 
shoulder like feature at k « 0.5257T, which is not too far 
off from the anticipated value of ir/2. In an ideally or- 
dered situation, it would be peaked like a delta function. 
But because the spin correlation decays rather rapidly, 
this shoulder like feature is all that we can see. Similar 
features for R = 6 occur around k = 7r/3 and 27r/3. For 
R = 8, they occur around k — tt/4, 7t/2 and 37r/4. So 
indeed, for a given even-i?, there are ^ — 1 non-trivial 
k points as anticipated from an Ising type argument. A 
better view of these non-trivial A-points can be gained 
by Fourier transforming (— l) r e~ r ^C{r) as shown in the 
insets of Fig. 6. It should also be noted that already for 
R = 8, the structure factor is almost like that of H 3 , as 
show in Fig. 9. 



FIG. 9. Comparing S(k) of H& with that of H 3 for L = 40. 
The two are almost identical. 
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FIG. 10. The nearest-neighbor concurrence vs. R. 

E. Entanglement 

We have also calculated some information-theoretic 
quantities such as the concurrence and the entropy (of 
entanglement) in the ground states of the even-i? mod- 
els. These are good measures of quantum correlation. 

Since the spins in the linear-exchange models are inter- 
acting via SU(2) symmetric antiferromagnetic exchange, 
the ground state in general is a singlet (total S — state) 
whether known exactly or not. For a given singlet state, 
the reduced density matrix of a pair of spins, obtained by 
tracing over the other spins, is also rotationally invariant, 
given by the following one parameter Werner state. 26 

P2(i,j)=p[i,j){i,j] + ^~I (4) 

Here, [i,j) is the singlet state formed by the spins Si and 
Sj, and p = — (4/3)(S^.Sj-). For p = 1, the Werner state 
is a pure singlet state. In general, however, it represents 
a mixed state. Therefore, the entanglement between two 
spins is given by the concurrence, C P2 = max(0, |p— \ ) . 2 ' 
For the even-i? linear exchange models, we have com- 
puted C p for the nearest-neighbor spins for L = 40. It 
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FIG. 11. Entanglement entropy plotted against the separa- 
tion r between the spins. The values of S(p2) at r = 1, which 
are very close to zero, are not shown here in order to zoom 
into the features at r ~ R that appear as small dips. 

is shown in Fig. 10, whose strikingly similarity to the 
ground state energy plot (in Fig. 1) is immediately noti- 
cable. It is not surprising however, because the nearest- 
neighbor pair state tend to approach the perfect singlet 
as even-i? increases. This we know from the dimer-dimer 
correlations. The C P2 simply follows this approach to- 
wards the singlet state, just as the ground state energy 
does. 

The entropy is another measure of correlation within 
a state. For a spin-1/2 pair, if the two spins are uncorre- 
cted, the entropy would be maximum, that is 2 (taking 
log base 2). In the Werner state, it is given by the rela- 
tion. 28 

S(p 2 ) = 2 - log 2 (l + 3p) 3^ log 2 (l -p) (5) 

Clearly, for a perfect singlet between a pair of spin-l/2s, 
that is p — 1, the entropy S(p2) is zero as it ought be. 
Figure 11 shows the S{p2) computed for different sepa- 
rations along the chain. Expectedly, it is close to zero 
for the nearest-neighbor spins. However, it quickly ap- 
proaches the value of 2 for larger separations, consistent 
with the rapidly decaying spin-spin correlations. Inter- 
estingly, while remaining mostly equal to 2, the entropy 
also shows signatures of slight recovery of correlations at 
a distance of i?, in accordance with the correlation-revival 
seen in the spin-spin correlation function. 

F. Evolution of energy gap as a function of R 

Consider a suitably modified version of the linear- 
exchange models, 

H B , a =H R + aV R , (6) 

where V R = J2i Sj. [Si+i + S; +2 + ... + Si +R ], and the 
parameter a varies from to 1. It has been noted 18 ' 19 



FIG. 12. Evolution of the spin-gap between successive linear- 
exchange models, as function of a defined in Eq. (6). 



that the successive i? models can be related through H R a 
as: H R+ i = H Rt i = H R + V R . The a can thus be viewed 
as a parameter which continuously changes H R to H R+1 . 
It offers an interesting handle to see how even-i? evolves 
to odd-i? and vice versa. Here, we present the evolution 
of spin-gap as a function of a (see Fig. 12). We find 
that the gap, A, remains non-zero (except for H2), and 
smoothly increases or decreases for a varying between 
and 1. The spin-gap for an odd-i? appears as a local 
maxima between two even-i? cases. 



IV. SUMMARY 

We have performed DMRG calculations on the linear- 
exchange quantum spin chains given by the Hamiltonian 
H R of Eq. (1). In particular, we have investigated the 
models for R = 4, 6 and 8, in order to understand the 
general nature of even-i? models in relation to the odd-i? 
cases (whose ground state properties have been known 
exactly). We have computed the ground state energy, 
the spin-gap, the dimer-dimer and spin-spin correlations, 
the nearest-neighbor concurrence, and the entropy con- 
tent in spin-pairs, to uncover the low-energy physics of 
the even-i? models. Based on our findings, we now confi- 
dently understand the class of linear-exchange models as 
summarized below. 

The integrable Hi stands out alone with algebraic spin- 
correlations in the ground state and gapless excitations. 
All the rest are found to be spin-gapped with rapidly 
decaying spin correlations. 

The ground states of the even-i? ^ 2 models exhibit 
nearest-neighbor singlet dimerization, akin to the odd-i? 
models. In fact, numerically the case of i? = 8 already 
shows the almost perfect MG dimerization. The rigorous 
asymptotic behavior, it appears, sets in rather quickly 
even for not-so-very-large values of even-i?. Thus, all the 
linear exchange models (except H 2 ) exhibit spontaneous 
dimerization like the MG model. 
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The one thing about even-i? ^ 2 models that is com- 
mon with H2 is the 7r-oscillations in the spin correla- 
tion (relics of the Neel type ordering in the correspond- 
ing Ising cases). While the spin-correlations decay very 
rapidly (for R ^ 2), they also show weak signs of (peri- 
odic) revival, before eventually decaying to zero. In the 
Fourier space, these features give rise to % — 1 identifiable 



fc-points roughly close to 2im/R (n — 1,2, .., § — 1), as 
suggested by the corresponding Ising picture. 
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